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Abstract: We investigate the scale on which 
the correlation arises between the 843 MHz ra- 
dio and the 60 fim far-infrared (FIR) emission 
from star forming regions in the Milky way. 
The correlation, which exists on the smallest 
scales investigated (down to ~ 4pc), becomes 
noticeably tight on fields of size 30', correspond- 
ing to physical scales of ~ 20 — 50 pc. The FIR to 
radio flux ratio on this scale is consistent with 
the radio emission being dominated by thermal 
emission. We also investigate the location de- 
pendence of gmoan, a parameter measuring the 
mean FIR to radio flux ratio, of a sample of star 
forming regions. We show that gmcan displays a 
modest dependence on galactic latitude. If this 
is interpreted as a dependence on the intensity 
of star formation activity, the result is consis- 
tent with studies of the Large Magellanic Cloud 
(LMC) and other nearby galaxies that show el- 
evated values for q in regions of enhanced star 
formation. 



1 Introduction 

The tight correlation between radio and far-infrared 
(FIR ) emission in exte rnal galaxies has been long known 
(e.g.. lYun et alJboOlh . To explain the correlation, the 
conventional picture suggests that both emissions are 
the result of massive star formation. The radio emis- 
sion has both a thermal component and a non-thermal 
synchrotron component due to cosmic ray electrons 
accelerated to relativistic speeds through supernovae 
shocks. The FIR is mostly thermal emission from dust 
excited by UV emission associated with star formation. 
The FIR emission is a common measure of star for- 
mation rates (SFRs, iKennicutt et al.|[r998l ), although 
a contributi on to the FIR e mission from older stellar 
populations l|Bell et alj|2003l ) necessitates a more com- 
plex calibration from luminosity to SFR than a simple 
linear relation. The effects of different dust geometries 
and opacities in the immediate vicinity of different star 
forming regions potentially adds a signific ant compli- 
cation (|Calzetti|[T997l : ICalzetti et al.l ll995l '). There is 
also a strong trend for higher levels of dust obscura- 
tion t o occur preferentially in ga l axies of higher SFR 
(e.g.. iPerez- Gonzalcz et al.l l2003l : lAfonso et all l2003l : 
iHopkins et al..,2001. ). so it is not unreasonable to ex- 



pect a similar variation in individual star forming re- 
gions. 

A complication in understanding the radio-FIR cor- 
relation arises from the fact that synchrotron radia- 
tion occurs on a larger spatial scale, due to its difi^u- 
sion length, than th e FIR thermal emission, (see, e.g., 
iMurphv etaD 120061 ). Measurements o f the total FIR 
to ra dio flux ratio, q llHelou et al.lll985l '). for the Milky 
Way l|Broadbent et alJll989l l show that on the largest 
scales our Galaxy is consistent with other star forming 
galaxies, although these authors also emphasise that 
the tight global correlation is not seen within galaxies. 
Measur ements of q within a few h undred pc around 
Orion (|Boulanger fc Peraultl Il988h give a value ele- 
vated compared to the average for star forming galax- 
ies, and similar to those seen in regions of active star 

for mation in the L arge Magellanic Cloud (LM C) by 

Hug hes et ahl ||2006D and within nearby galaxies (jMurphv et al.l 
,2006 ). A question of interest then becomes, on what 
spatial scale does the radio-FIR correlation begin to 
appear around star forming regions? It has been pro- 
posed that this should start to appear on scales compa- 
rable to t he Ri 1 — 2 kpc d iffusion length for cosmic ray 
electrons (|Bicav fc Helou 1990: Marsh. & Helou 1998). 
This would correspond to the observed region around 
an area of massive star formation being large enough to 
include sufficient quantities of radio synchrotron emis- 
sion. We show that, for the Milky Way, a radio-FIR 
correlation exists on much smaller scales, and is likely 
to be a correlation between the thermal radio and FIR 
components. 

Recent studies of the radio-FIR correlation have 
concentrated on analysing the physical processes which 
lead to the correlation on galactic scales. This has 
been d one by analysing the correlation on small spatial 



scales l|Xu et al.ll 19921 : iHughes et al 



other nearby galaxies (jMurphv et al 



2006 ). and within 



20061 ). iHughes et all 



(|2006l ). looking at the Large Magellanic Cloud (LMC), 
achieved a spatial resolution of 20 pc. On such fine 
spatial scales, these authors have demonstrated vary- 
ing q for different structures in the LMC. In regions 
of the most active star formation, q is higher. This 
is interpreted as due to two different physical mecha- 
nisms driving the radio-FIR correlation. In the most 
active star forming regions, radio emission is mostly 
thermal while in less active regions, the correlation is 
due to synchrotron radiation. We compute the mean g, 
called Qmean herein, for star forming regions grouped by 
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their galactic latitude. We show that gmcan is elevated 
closer to the galactic plane, where the most active star 
formation in our galaxy is occuring, compared to the 
value at higher galactic latitudes. 

2 Analysis 

We investigate the spatial scale on which the radio-FIR 
correlation begins to arise around massive star form- 
ing regions in the Milky Way. We then extend this 
analysis to explore whether the radio-FIR correlation 
shows any Galactic latitude dependence, as might be 
expected if q depends on star formation rate, as sug- 
gested by th e elevated q values seen in the LMC bar 
( Hughes et a l. 2006) and spiral arms of nearby galaxies 
( Murphv et al.ii200Ql . 

2.1 Data 

The star forming regions used for this analysis are 
drawn from the Census of High and Medium-mass Pro- 
tostars (CHaMP) survey conduct ed with the Mopra 
dish o f the Australia Telescop^H (^Barn es et al ] l2006l . 
l2009al '). CHaMP targeted m assive dense gas clum ps 
based on the Nanten Survey l|Yonekura et al.ll2005r ) of 
carbon monoxide gas, a standard tracer of the sites 
of star formation. The CHaMP catalogue is based on 
an examination of complete Nanten maps for the re- 
gion 280° < I < 300°, -4° < fe < +2°, made in four 
J = 1 — molecular spectral lines in the 3 mm band: 
CO, / ^CO, C^^O, and HCO+ (see also lYonekur"a et al.l 
l2005l . for a subset of these data). The CHaMP cata- 
logue contains 209 local emission maxima in {l,b,v) 
space, ident ified from the Nante n C'^^O and HCO"'' 
data cubes l|Barnes et al.l [2009bl ). In this sense, the 
CHaMP catalogue comprises an unbiased sample of all 
the dense gas clumps in this window, and is therefore 
a useful indicator of the most likely sites of current 
massive star formation in this part of the Milky Way. 
The unbiased nature of the CHaMP selection, allow- 
ing identification of all massive star forming regions 
throughout a representative spiral arm of the Milky 
Way, suggests that our analysis should be generally 
applicable to the Galaxy as a whole. 

Using the coordinates of the 173 highest CO lumi- 
nosity sources from the CHaMP survey, being those 
already catalogued at the initiation of this study, 138 
target regions were selected for this analysis based on 
availability of both radio and FIR images. The omis- 
sion of the 20% of targets lacking radio and FIR data 
should not introduce any significant bias to our analy- 
sis, as we still sample the full range of Galactic latitude 
and longitude within the CHaMP survey area. The ra- 
dio images at 843 MHz were obtained from the Sydney 
Univ e rsity Molonglo Sky Survey (SUMSS. lBock et all 
Il999l : iMauch et all l2003l 1 and FIR images at 60/im 

^The Mopra telescope is part of the Australia Telescope 
which is funded by the Commonwealth of Australia for op- 
eration as a National Facility managed by CSIRO. The Uni- 
versity of New South Wales Digital Filter Bank used for the 
observations with the Mopra Telescope was provided with 
support from the Australian Research Council. 



from the Infrared Astronomical Satellite (IRAS) Sky 
Survey Atlas (ISSAfl An advantage of the SUMSS ra- 
dio images for this analysis is the high sensitivity to 
extended emission, a consequence of the closely-spaced 
and highly redundant spacings between rec eivers in the 
Molonglo Observatory Synthesis Telescope (jBock et al.l 
1 19991 ). The resolution of the SUMSS images is 45" x 
45"cosec|(5|, and the internal flux densi ty scale is ac- 
curate to within 3% (|Mauch et al.l[200^ . An example 
of both SUMSS and IRAS images is given in Figure [TJ 
showing a region two degrees on a side, centred on one 
of the CHaMP targets. 

2.2 The spatial scale on which the 
correlation arises 

For each CHaMP star forming region, SUMSS and 
ISSA images of varying sizes were obtained. Square 
field sizes of 5', 10', 15', 30', and 1° were analysed. Re- 
gions 2° on a side were also obtained, but most of the 
SUMSS images of this size were incomplete. This is a 
limitation of the SUMSS image server, which does not 
mosaic adjacent images if the requested field extends 
beyond the boundary of the nearest 4° x 4° archival 
image. To measure the correlation between the FIR 
and radio images of the same star forming region, a 
simple flux sum was taken by adding up the ampli- 
tudes of all pixels within the image. Figure [2] shows 
radio fiux density as a function of FIR flux for the four 
smallest spatial scales measured. Each of the 138 star 
forming regions is shown, except for those flelds where 
the fiux summed to a negative value in the radio im- 
ages due to the presence of artifacts. Such fields were 
not included in any of our analyses. 

The correlation coefficient was calculated for each 
field size (Table [I}. A significant correlation is found 
for all, even the smallest field sizes. By fields of size 
30' X 30' this correlation becomes visibly tight, espe- 
cially for the brightest sources (associated with high 
SFRs). The correlation at the smallest scales (5' to 
30') was tested using FIR images randomly offset from 
the radio images. The correlation coefficients found for 
these non-associated regions are significantly reduced, 
as would be expected if the correlation is due to the 
emission from the massive star forming regions. For 
the randomly offset fields, the correlation coefficient 
varies from 0.61 for the 5' fields to 0.64 at 15', increas- 
ing slightly to 0.75 at 30'. This supports the suggestion 
that, even on the smallest scales measured here, there 
is a measurable radio-FIR correlation associated with 
the CHaMP massive star forming regions. A value 
for Qmcan was Calculated, along with the corresponding 
rms scatter in the individual q values for each field. 
The observed rms scatter in q for our targets is con- 
sistent with t he ranges of q foun d withi n gala xies by 
iHughes et al] (2006) and Murph v et al] (|2006l ). It is 
likely that much of this observed scatter is real and 
not simply a result of limitations in our approach, or 

2 This research has made use of the NASA/ IPAC In- 
frared Science Archive, which is operated by the Jet Propul- 
sion Laboratory, California Institute of Technology, under 
contract with the National Aeronautics and Space Admin- 
istration. 
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Figure 1: SUMSS 843 MHz (left) and IRAS 60 (right) images for a 2° x 2° region centred on one of 
the target CHaMP massive star forming regions. 
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5 


0.86 


3.41 


0.54 


10 


0.83 


3.34 


0.54 


15 


0.87 


3.33 


0.47 


30 


0.87 


3.42 


0.48 


60 


0.81 


3.56 


0.44 



Table 1: The correlation between the FIR and 
radio emission around regions identified in the 
CHaMP survey in the Milky Way. Values of gmoan 
and the rms scatter in individual q values are also 
shown. 



due to measurement uncertainty. The carbon monox- 
ide isotopologues, used as a tracer by tlie CHaMP sur- 
vey, identify the locations of the coldest and densest 
gas clumps. Such dense gas clumps are variously asso- 
ciated with pre-stellar gas, active star forming regions, 
and residual clouds from recent massive star formation. 
Such variation could account for a significant portion 
of the scatter observed in q. 

In order to ensure these measurements provide ro- 
bust estimates of Qmcan, we explored several possible 
sources of bias: (1) The presence of artifacts in the 
radio images; (2) The contribution of a non-zero back- 
ground level in the ISSA images; (3) Overlapping fields, 
arising from the clustering of CHaMP sources. For 
each of these tests the values of q and Qmean were re- 
measured using the 30' x 30' images, and the corre- 
lation coefficient recalculated to ensure no bias was 
introduced. The initial correlation coefficient found 
for the 30' x 30' images was 0.87. Images removed 
in the artifact analysis include those displaying ra- 
dial artifacts originating from bright sources, and those 



with large negative-flux regions associated with bright 
sources, along with images that did not encompass the 
full field of view as a consequence of the SUMSS im- 
age archive limitations mentioned above. Excluding 
those radio images visually identified as containing ar- 
tifacts marginally decreases the correlation coefficient 
to 0.84. It was noted that the rms scatter of q in 
the sample with visually identified artifacts removed 
was not significantly changed, indicating that such ar- 
tifacts are unlikely to be the primary source of the ob- 
served scatter. The median background level in each 
IRAS image was measured from the histogram of pixel 
intensities, and subtracted prior to calculating gmcan- 
This resulted in only minor changes to the resulting q 
values as the background level was low compared to 
the total flux. The resulting correlation coefficient is 
0.85. To identify the impact of image overlap, (/mean 
was calculated only for images independent by more 
than 75% of their size. The resulting correlation co- 
efficient is 0.74, possibly suggesting that overlapping 
fields account for some of the strong correlation seen 
in Table[T] A smaller fraction of the 15' fields, though, 
are similarly overlapping, and removing those does not 
significantly change the correlation coefficient found 
(0.86) compared to the full sample (0.87) for that field 
size. We conclude that our estimate of the field size on 
which the radio-FIR correlation arises is not sensitive 
to these limitations. 

To infer the spatial scale associated with this an- 
gular size we require the distances to the CHaMP star 
forming regions. Although there are currently only a 
limited number of distances available, we make use of 
those in hand to estimate the length scale on which the 
radio-FIR correlation arises. For our sample of 138 
star forming regions, there are 54 distance measure- 
ments, based on a combination of techniques. Some of 
the objects have distance estimates based on nearby or 
embedded HII regions, while most are estimated kine- 
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Figure 2: SUMSS 843 MHz radio flux density as a function of IRAS 60 /im flux, for regions of progressively 
larger size centred on CHaMP massive star forming regions. The radio-FIR correlation tightens noticeably 
once the area analysed reaches about 30' x 30'. This corresponds to a physical scale of « 20 pc to 50 pc. 



matically. The {l,v) diagrams in the direction of our 
sample fields are largely unambiguous, with distance 
uncertainties of ~ 10 — 20%. Less than 50% of the dis- 
tance measurements suffer from a near/far ambiguity. 
Using these estimates, we can identify the approximate 
range on which the radio-FIR correlation begins. An 
angular scale of 30' corresponds to a spatial scale of 
22 pc to 59 pc at the estimated distances of 2.5 kpc to 
6.8 kpc. This reflects an approximate diameter cen- 
tred on the star forming regions. The 5' field sizes, 
on which a correlation also exists, correspond to (di- 
ameter) scales of 3.6 pc to 9.9 pc. These scales are all 
significantly smaller than the dif fusion length scale ex - 
pected for cosmic ray electrons (|Bicav fc Heloul [19901 ) ■ 
and support the idea that correlations between differ- 
ent components of the radio and FIR emission exist 
over a variety of spatial scales. This result will be ex- 
plored further in the discussion. In the next section 
we quantify our measurement of gmcan and explore its 
galactic latitude dependence. 



2.3 Galactic latitude dependence 

The average 60/im FIR to 843 MHz radio ratio of se- 
lected images is measured by the ijmcan parameter and 



is defined as: 

gmcan = (logio ( -^^^^^^ )>• (1) 
V <5843MHz / 

Note that this differs from the q typically defined in 
that here we use different FIR wavelength and ra- 
dio frequency. This definition of gmcan will not af- 
fect the trends observed. The use of 843 MHz rather 
than 1.4 GHz radio measurements will act to reduce 
our gmoan by up to 0.15 ( for pure synchrotron emis- 
sion) compared to that of iHughes et al.l l|2006D . This 
estimate simply assumes a conversion from 843 MHz 
to 1.4 GHz using a spectral index a = —0.7 (with 
S cx v"). Using the 60/im flux only, rather than the 
combined 60 /xm and 100 /xm FIR flux reduces our esti- 
mate of (/moan by approximately 0.15 again, estimated 
assuming that on average 5 'ioou ~ 1.6S' 6 0mn, a nd us- 
ing the definition of FIR of iHelou etal] (| 19851 '). This 
total systematic offset of « 0.3 is slightly smaller than 
the rms scatter in our measurement (Table [T]). 

Values for Qmean are shown in Table [T] After cor- 
recting the FIR images by subtracting the background 
level, we find gmcan ~ 3.34 for the 30' fields, and 
Qmcan = 3.50 for the 1° fields, only marginally smaller 
than the initial estimates, given the observed scatter. 
The values that we find here for (/mean, especially af- 
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ter considering the systematic offset of ~ 0.3 indicated 
above, correspond to flux ratios more tlian an order 
of magnitude larger tlian those corresponding to the 
typic al value of qr ~ 2.3 found for star forming galaxies 
fe.g.. ICondonllTOQ^ ). Our estimates are, however, con- 
sistent with the elevated values of q identified within 
nearby galaxies in regions of elevated star formation 
dHughcs ct al. 2006; Murphv et al. 2006). 

To measure any possible galactic latitude depen- 
dence of Qmean, we use the 1° images, although the re- 
sults are similar with other field sizes. The images were 
divided into two bins of galactic latitude, \b\ < 0.8° 
and \b\ > 0.8°, containing 39 and 32 target regions re- 
spectively. The value of gmoan shows a clear galactic 
latitude dependence, having higher values on average 
closer to the galactic plane (Figure A Kolmogorov- 
Smirnov (KS) test comparing the gmoan values for the 
two subsamples gives P < 1 x 10^'^, indicating that 
they are inconsistent with having been drawn from 
the same population. The Seo/jm fiux values for each 
of these subsamples were also compared using a KS 
test. This shows again that, with P < 1 x 10"'', the 
60 /xm fluxes of the subsamples are inconsistent with 
being drawn from the same population, and that those 
at lower latitudes show higher 60 fim fluxes consistent 
with having higher star formation rates. 
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Figure 3: Values of g as a function of S'eo/^m for the 
two subsamples split by distance from the Galac- 
tic plane, gmoan is elevated closer to the Galac- 
tic plane, where the star formation, indicated by 
SeOfj.m, is most active. 



3 Discussion 

The first main result of this investigation is that a local 
radio-FIR correlation seems to arise at spatial scales of 
at least « 20 — 60 pc, and likely on scales as small as a 
few pc. The second main result is that the FIR-radio 
ratio (gmcan) has a galactic latitude dependence. We 
interpret this second result to indicate q is dependent 
on the level of star formation activity. 

The scale on which the radio-FIR correlation ap- 
pears around star formation regions is consistent with 



the LMC study (|Hughes et aLlbOOeT ) where the correla- 
tion arises at ~ 20 pc and becomes tight above ~ 50 pc. 
We further find a correlation on even smaller scales, 
down to ~ 4pc. The diffusion length of relativistic 
electrons co ntributing to th e galactic FIR-radio ratio 
(« 1 - 2 kpc, iBicav fc Heioul fl990') is much larger than 
these spatial scales. This suggests that the correlation 
on such a local scale is between thermal radio and FIR, 
an interpretation consistent with other studies show- 
ing strong corr elation between thermal radio and warm 
dust emission (|Hoernes et al .ll998l : lHugh es et al.ll2006l l. 
The the rmal frac t ion at 843 MHz calculated from Equa- 
tion 5 of lCondonI (|1992| ') is about 0.08, so assuming that 
the FIR emission is constrained to the locale of the star 
forming region, the fiux ratios observed would be ex- 
pected to be about an order of magnitude larger than 
those corresponding to g ~ 2.3 found as the average 
for star forming galaxies in general. This is consistent 
with the values that we have measured, supporting the 
result that the high value s of (7mean from this study, a s 
well as earlier estimates teoula nger fc Peraultl ll988''). 
are a consequence of not including the synchrotron ra- 
dio emission that is distributed on larger scales. 

The result of the latitude dependence analysis is 
consistent with a dependence of gmean on the star for- 
mation rate. This is based on the assumption that star 
formation rates are higher closer to the galactic plane, 
supported by the increased Seo/jm measured there. Our 
results are consistent with those o f lHugh es et al.l (|2006t ). 
who observed elevated q values in the bar of the LMC 
where there is elevated star formation activity. Strong 
additional support for our intepretation of the galac- 
tic latitud e analysis comes from a study of local spiral 
galaxies. iMurphv et al.l l|2006l l observed that active 
star forming regions of galeixies with high SFR display 
elevated q values. 

From Table [1] the change in gmcan between the dif- 
ferent spatial scales is unlikely to be statistically signif- 
icant given the rms scatter in the measurements. Even 
out to scales of 2° (spatial scales of up to ~ 200 pc), 
where these were able to be measured, the value of 
gmcan remains essentially unchanged, although there 
is perhaps a very marginal elevation on scales of 30' 
to 1°. This, while not statistically significant, may 
imply a suggestive trend. We suggest that a rise in 
gmean at « 20 — 50 pc aud the decrease on larger scales 
around star forming regions could be due to the vari- 
ous components of FIR and radio emission occuring on 
different spatial scales. The change in gmean as spatial 
scale increases is then a consequence of the changing 
proportions of these different components. We can de- 
scribe such a behaviour with a simple model for the 
distribution of these components consistent with the 
known spatial scales of each emission process. 

On the smallest scale, thermal FIR and radio are 
dominant. The thermal radio emission is due to free- 
free emission in ionised gas while the thermal FIR 
emission is from hot dust. On a slightly larger scale, 
cirrus FIR emission continues to contribute, eventu- 
ally ceasing and leaving only the synchrotron radiation 
component on the largest scale. 

The simple model shown in Figure [4] demonstrates 
that any rise in gmean between the 30' and the 1° im- 
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Figure 4: A simple model of the radial scale on 
which the components of FIR and radio emission 
extend, centred on a star forming region. 



ages could be due to the cirrus FIR emission contin- 
uing to greater spatial scales than the thermal radio 
emission. The value of q, of course, is calculated using 
fluxes corresponding to the integral of these intensities 
within a given scale. On the kiloparsec scale, and for a 
single star forming region, q would continue to decrease 
with the increasing contribution of the synchrotron ra- 
diation, since all the FIR emission has ceased. Within 
galaxies, the overlapping of adjacent star forming re- 
gions on scales comparable to the synchrotron diffusion 
length then gives rise to the observed lower average 
value of q dominated by the synchrotron component. 

This model, focussing on the dependence of q on 
spatial scale, does not directly address whether or not 
q should be higher around a region of higher star for- 
mation rate, for the same spatial scale sampled (Fig- 
ure [S]). This may be the case, for example, if the rel- 
ative proportion of thermal FIR and radio emission is 
itself dependent on the star formation rate. 

4 Conclusion 

We have examined the scale on which the radio-FIR 
correlation arises around star forming regions in the 
Milky Way. The correlation is tight by scales of ~ 20 — 
60 pc, and exists on scales even as small as « 4pc. 
The galactic latitude dependence of gmoan is consistent 
with that found in the LMC and local spiral galaxies 
(|Hughes et al.ll2006l : lMurphv et al.ll200^ ). and strongly 
suggests that q is elevated in more active star forming 
regions. We interpret this as a correlation between the 
thermal radio emission and the FIR. 

We have suggested a simple model that accounts 
for these observed trends in gmcan- To test this model, 
we plan to compare individual q values, on a range of 
spatial scales, with the directly measured star forma- 
tion rates in each CHaMP target region. Star forma- 
tion rates can be measured in a number of different 
ways, including the radio and FIR luminosities mea- 
sured specifically to encompass each star formation re- 
gion. For such radio measurements, which are likely 



to be primarily thermal emission, existing star forma- 
tion rate calibrations for total radio luminosity can 
be revised b y referring to the thermal to total radio 
flux ratios of lCondonl (|l992h . Ha luminosities, as well, 
may be able to be obtained for many of the regions, 
and used in estimating star formation rates. Addi- 
tional data for the CHaMP targets is also planned to 
be obtained through programs on Gemini-S, SOFIA, 
and Herschel. 

Multifrequency radio continuum measurements for 
each star forming region would allow determination of 
the thermal-to-synchrotron ratio as a function of spa- 
tial scale. This would provide a more accurate quanti- 
tative constraint on the scales within the simple model 
proposed here. It would also allow a refinement of the 
radio-derived star formation rates. 

We warmly thank the referee for extensive and 
thorough comments that have led to significant im- 
provements in this paper. AMH acknowledges support 
provided by the Australian Research Council through 
a QEII Fellowship (DP0557850). PJB acknowledges 
support through NSF grant AST-0645412 at the Uni- 
versity of Florida. 
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